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Мета статті - отримати математичні конструкції для геологічних об'єктів типу синкліналей та антикліналей, об- 
грунтувати єдиність оберненої задачі відновлення аналітичних моделей горизонтально-шаруватого геологічного середо- 
вища з кількома густинними межами розділу для цих конструкцій у визначеному заздалегідь класі Чорного контактних по- 
верхонь та апробувати розроблену методику для їх ітераційного обчислення. Сукупність цих двох моделей утворює нову 
уточнену постановку оберненої задачі гравіметрії для контактної поверхні. Це необхідно для покращення відомих процедур 
підбору у розв'язанні обернених задач гравітаційних і магнітних полів. 

Обернена задача визначення контакту у горизонтально-шаруватому середовищі з кількома густинними межами зведе- 
на до розв'язання нелінійного інтегрального рівняння, яке описує контакт, обмежений заданими сталими асимптотами у 
плоскій області. Проте у такій постановці практика обчислень ускладнюється проблемою еквівалентності розв'язків. 

Для цієї моделі наведені дві теореми розділення полів - для випадку кількох однозв'язних об'ємів і для випадку кількох 
шарів, що не перетинаються. Теореми єдиності засновані на теоремах розділення полів, які дозволяють звести розв'язок 
оберненої задачі за сумарним зовнішнім полем п об'єктів (рудних тіл, меж розділу шарів) до розв'язку оберненої задачі для 
окремих об'єктів - за значеннями поля від цих геологічних об'єктів. 

Вказано чисельні схеми для визначення початкового наближення густинного контакту у багатошаровому геологічному 
середовищі. Ці алгоритми формально співпадають на першому кроці ітерацій. Аналогічні схеми на основі ітераційної конс- 
трукції Чебишева запропоновані і для ітераційного уточнення поведінки "асимптот" контакту. 

Здійснено моделювання синтезованих початкових наближень антикліналей" і "синкліналей" за цими алгоритмами. Вка- 
зано альтернативний спосіб обчислень, який базується на визначенні різних моментів кривої контакту. Для обчислення 
інтегралу отримано відповідний вираз у скінченних квадратурах. 

За результатами моделювання виявлено, що нові аналітичні конструкції для обчислення багатошарових контактів 
при їх чисельному моделюванні за способом Ньютона швидше збігаються порівняно із класичними методами обчислень 
контакт:у. Їхню стійкість на даних великої розмірності доцільно перевірити на польових даних. Спроби обійтися грубими 
наближеннями успіху не мали: збіжність на порядок менша та досить сумнівний геологічний зміст. 

Ключові слова: теорія потенціалу, аналітична модель, контактна задача, класи контактних поверхонь, розділення гра- 
вітаційних полів, моделювання. 


Постановка проблеми. У рамках перебудови ана- свячена ця стаття. 
літичної бази теорії потенціальних полів визначено Загалом єдиність розв'язку обернених задач граві- 
аналітичні моделі гравітаційного поля та складно по- метрії гарантована лише в кількох класах аномальних 
будованого  горизонтально-шаруватого геологічного джерел. Зокрема, для контактних задач однозначне 
середовища |5). Перша модель отримана з еквівален- визначення густини 5(5и ба) і кусково-неперервної межі 


тного зображення рівняння сили тяжіння |/), який ви- 


іл вищ х праведливе для кл т- 
пливає з аналізу плоского нормального потенціалу розділу середовищ х, (5,) справедливе д асу Ос 





сили тяжіння у локальному околі точки вимірювань. ромогильського Оз" (а,б) ,у якому 

Хоча можливий і інший підхід до розробки аналітичних 

конструкцій поля |9). Друга модель формалізована х(е)- -п- бх, (5,), 5 є|а,Б зх, (а) - бх. (6)- 0 (1) 
обмеженнями, властивими класу Чорного контактних щі -В, сє|абі" з ; 
поверхонь. Ці моделі виросли з нових підходів до чи- РЕ 

сельного обгрунтування методів обробки великих об- зх, (х | 2с(х, -а) ,асххаєб а» 2 
сягів геофізичної інформації. Сукупність цих двох мо- б (х) а оку, вбавуявіві зад 

делей утворює нову постановку оберненої задачі 

гравіметрії для контактної поверхні. та Страхова 54, Б), де 


Аналіз публікацій. Визначення контакту у відомих 
моделях середовища |8| трактується як обчислення хів)з| п я дж (5,), 5, є Га) | 5х (а)- 0, йо 20, (3) 
різниці відхилень деякої теоретичної моделі середо- р 





вища від заданого її опорного елемента, від якого ці ів (Є) за(5 -а) і асржазб азо 
відхилення не надто великі. Але при цьому постулю- (4) 
й й фах паз ЛИ , 

ється наявність апріорної інформації про елементи зх, (а 2 6 (-5,) й ребивьо бо 
геометрії (глибину, форму, орієнтацію) цього опорного ти пр 
елемента моделі для формування початкових набли- сі, с, Є|а, | - внутрішні або граничні точки відрізку 
жень. Вдалий вибір початкових наближень досі визна- (а, Ь1, ана нескінченному профілі 
чає успіх ітераційних алгоритмів підбору |2). М ана 

Для горизонтально-шаруватого середовища з кіль- 7 - (5) -ПОКО(С), 5 Є (-бо» бо), 
кома густинними межами шарів введено аналітичні мо- 2(2)- 1, « а, ЇЇ", ІЗ 2 (5) 
делі поля і геологічного середовища, в яких опорні | 

2(2)-р «561, в»о0 

елементи (глибина і асимптоти) розраховуються в про- 2(5)- М, «6,5, 8»0, 
цесі ітерацій. Ці розв'язки у вигляді нелінійних рівнянь 
справедливі ,в малому", для локальних областей, роз- де б, В,2»0. Тут 28 7 2(с,) є 8Ир хі 2(5,) | 
міри яких істотно залежать від розмірів аномальної об- 
ласті і точності апроксимації складових поля |41. й 2. - 2(с,)- п! ехії (69 при с с, є | -5., 5, 

Мета статті. Обгрунтуванню єдності цих моделей у Ї-8о бо 
визначеному заздалегідь класі СП" (-со, с) та апро- -2, «Р, -2, 40, б« ПВ, «Й, 


бації конструктивних алгоритмів для їх обчислення при- 
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Зауваження 1. Підкласи Ов і Ов класу 


Озі9 (аб) можна вважати моделями антикліналей | 


синкліналей. 

Загалом задача визначення контакту у горизонтально- 
шаруватому середовищі з кількома густинними межами 
зведена до розв'язання нелінійного інтегрального рівнян- 
ня, яке описує контакт, обмежений заданими сталими 


асимптотами у плоскій області С, є СП'"Р(-о0, со) |БІ: 


р 2 2 
-хузчЗН 
шодо кіп КН аа, (6) 
а (5, -ХИ) - Х5(5,) 
Однозначність відновлення контакту у цій задачі га- 
рантована на класі Чорного СП"? (-оо, с) за умов, що 


значення ЛП (скінченна дуга кривої 2(2)) і величина 
стрибка щільності 25 - 6, - б, відомі. Клас СП ( о, со) 
визначено, як підклас Страхова, а саме, у класі (5) ви- 
ділена горизонтальна смуга С , обмежена згори і знизу 
прямими 2, «й, 2-й,, 0 «й, «Й, «чоо. Цю смугу роз- 
діляє на два шари із густиною 6,, 5, нескінченно глад- 


ка крива СС" з асимптотою 2-й, П,«П«й,, зірчас- 
тою щодо нескінченно віддаленої точки. На смугу дода- 


тково 2(2)-8|«КІЇ,  а»б, 


З 3 бу г(б)є дб. 


Така геометрія множини коректності дозволяє без 
зайвих ускладнень формувати початкові наближення 
для обчислення множин контактів. Проте у такій поста- 
новці практика обчислень ускладнюється проблемою 
еквівалентності розв'язків. Проілюструємо цю проблему. 

Проблема розділення поля. У практиці гравірозві- 
дки вирішальну роль відіграють задачі, у яких зовнішнє 
поле створюють складні розподіли мас, наприклад, п 


попарно неперетнутих об'ємів С,, К-1,2,..,п, або п 


накладена умова 


меж розділу шарів різної густини, що попарно не пере- 
тинаються. У задачах із складними розподілами мас 
єдиність розв'язку оберненої задачі вірна для досить 
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вузьких класів, тобто, за наявності великої апріорної 
інформації про джерела аномалій. 

Теореми єдиності оберненої задачі гравірозвідки 
для складних розподілів мас засновані на теоремах 
розділення полів, які дозволяють звести розв'язок обе- 
рненої задачі за сумарним зовнішнім полем п об'єктів 
(рудних тіл, меж розділу шарів) до розв'язку оберненої 
задачі для окремих об'єктів - за полями від цих об'єктів. 

Приклад 1. Нехай зовнішнє поле створюється ма- 
сами з носієм у вигляді об'єднання п скінченних однозв'яз- 


них об'ємів СЄ,, К-12,..,п, що попарно не перетина- 
ються, з межами О0С,, що не мають попарно загальних 


точок. Густина мас в кожному з об'ємів С, - обмежена. 

Теорема розділення полів. Якщо відомі зовнішні 
один щодо одного гладкі замкнуті кінцеві поверхні 5,, 
К-142,...п, строго усередині кожної з яких є одна і 
тільки одна компонента б, носія мас, то за заданим 
сумарним зовнішнім полем (поза об'єднанням об'ємів 
С, ) поля компонент б, визначаються однозначно. 


З цієї теореми легко отримати теореми для склад- 
них розподілів мас, аналогічні теоремам Новікова і Ра- 
ппопорта-Сретенського -- для об'ємів С, , зірчастих що- 
до заданої точки, чи таких, які мають середню густину. 
Густина ос(2) і система контактних меж, що не перети- 


наються, відновлюється однозначно |ЄЇ. 
Приклад 2. Нехай в п-шаруватому двомірному 


середовищі кожен з шарів С, має постійну густину 9, . 
Нехай межі розділу шарів 26, - нескінченні криві з го- 
ризонтальними асимптотами 


2-й і 0 сей, 


(Ак 20, В. » 0) ліворуч і праворуч (рис. 1). 





С сей 
Пі зба 1 
т є 
2 а б ретжужвоможно в анолн десни оце оре аура вва 
й- о б - -черііь "рибі рів ор бін У я тро ер чор, р, тя ро БАН Р у у и З б ч- 

чу 


ХЗеН 


Рис. 1. Розподіл густин у декількох шаруватих джерелах аномалій 


Теорема розділення полів. Якщо кожна із меж 
розділу С, містить горизонтальний відрізок Об, 


(5 - соп5і , а, 2 5 2 б, ), відомі проекція деякої його час- 
тини на вісь Ох а«, «Х «В, і прямі у,, які відділяють 
межі шарів, то за заданим зовнішнім (при 2 «0 ) сума- 
рним полем Лд з точністю до сталих визначаються 
поля 29, від окремих меж. 

Наслідок |3). Якщо при цьому відомі положення 
асимптот 2-й і 2-й; і перепад густини б6,, межі 


с, (К.1,2,... п) розділу шарів відновлюються за 


заданим сумарним полем однозначно. 

Задля відновлення контакту на основі виразу (б) за- 
пропоновано ряд ітераційних схем (16)-(18) 15). Розв'яз- 
ки легко отримати за відомого 2,(2): 2,(5) - 2,(5)-- В, (5). 
Зауважимо, що ці алгоритми формально співпадають 
на першому кроці ітерацій, оскільки 2, - 2, а на насту- 
пних кроках вони - формально різні. Аналогічні схеми 
запропоновані і для ітераційного уточнення поведінки 
"асимптот" Н-осоп5і. Власне, початкове наближення 


контакту можна задають як 2,(2) 2 -Н., - 49(5,Н,) або 
2(5)2-Н, ЧК49(5), деасосср, Н,»0 - задана ста- 
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ла ліва "асимптота" функції 2, Є |а, Б|, 


49(5) - 9(2)-9(а), 49(5,Н,) - аналітично продовжене 
на рівні -Н, значення Д4(25), К, - шукана амплітуда 
ундуляцій. 
Для обчислення К, визначена схема Чебишева |5 
" 2 

Р Ю ІНВ СО ОО (7) 
Й(2,) 2" (2) 
де 

2) « (1-сут((1-с)у 27) 


7 
-с п(с" - 27)--22 агсід-5------- 
7 реч 


с (1-с) 
-(1-суп(І- су -«Н71- с пс" «Н?)- 
-2н ага ос 9. (8) 


Зауважимо, що послідовні наближення кореня 7"? 


рівняння  7(2)-0 можна дотримати простіше, ніж 
П. Л. Чебишев. А саме: якщо корінь 2"? розташована у 
вузькому інтервалі |а, ВІ, і 2, обрано звідси ж, то з то- 
чністю до малих 3-го порядку щодо 708-2,; 77), 
2, Є| а, В| маємо: 

Р(2,) 

2 

Звідси при (2, )- 0 (що вірно, бо лише для кратних 


0 з (219) а (2) (2, 29-22) (219 - 2. У. 





коренів 7(2) - (2) - 0 ) отримуємо 
(2,)  (2,) 














2100 2 -- (29 -2 У. 
й Р(2.) 2Р(2,) й 
Завдяки | малості різниці | (2-2), поклавши 
2 " 2 
(20) -2) зр 2.) , отримаємо 70) -2 зе Н2,) Я (2,)ї 2) 
Р(2,) (2) 2Р(2,)) 


звідки випливає вираз (7). 
Зауваження 2. Величини Г(2.)20, М(2)20, а за 


значних Но» 0 маємо (2.)2 0, Г(2)«0. Отже, корені 
72в, Ху рівнянь (2,,9-)-0,(2,,9)-0 дають набли- 
ження 2. і 2, знизу і зверху, що впливає на "якість" 
початкового наближення 7, - 2, (5). 

Зауваження 3. Величини 2. і 2, невпевнено визна- 


чаються при великих Н , особливо якщо самі 7. і 2 





аг МА 1 АУРА, й 
" - я 
| .7 - п 


МИ ї 
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руувненнам а Зах  ЛРОРІНААНАРАРЛЕСНЕНИ ЕОЕРЕРИВИ НИ НРИНЕ ЗЕРЕН ЕРНРЕНИР ОВ 


- "То о 


великі. Тому додатково, крім 2. і 2, визначають при 
99 - 9(Б)-4(а) - 0, крім заданої лівої Н і праву "аси- 
мптоту" Н" кривої 2 - 2(2), тільки у виразі Г(2) замість 
9 або 9, беруть 59 для визначення Н". 


Враховуючи всі припущення попередньої статті |5| 
та зроблені вище, приступаємо до обчислень. 

Моделювання початкових наближень. Розглянь- 
мо схему відшукання початкового наближення для кон- 
тактів, що описуються функціями (3, 4) з класів "анти- 


чн 


кліналей" Оз, (а, б) і "синкліналей" Оз (а, Б). Ти- 
повими представниками цих класів є функції типу ,ан- 
тикліналь" та ,синкліналь", відповідно (рис. 2): 
22)--Н-а(1-2)8, 0421, (9) 
2(2)--Н-а(1-2)85, 0х«2 1. (10) 
Для відшукування К, систему координат зручно обрати 
так, щоб точка екстремуму д(х) була на початку коорди- 
нат, тобто ехіг 4(2) - 9(0) . Тоді і інтервал |с, 4) обираєть- 
ся "симетрично" щодо початку, так само як і а--с, р-с. 
При цьому для моделі (9) (ліворуч на рис. 2) вірні відношен- 
ня а, з- агсід С/Н, На, «га,, а» а,, а для моделі (10) 


а, - агсідС/2, На, » 20, а, «ОЦ». 
Функція (2) - 2сіп(с" - 77)- 42 агсід- -2с Іп(с" --Н?)- 
7 


С МДР 
-4Н агсід 19. де 9-9: для антикліналі, 9-9; - 


для синкліналі (обидва значення - зі своїми знаками). Її 








похідні Р(2)- 4агсід 3; Р(2)2- ра: Й Оскільки 
7 с" 
"(Н)-д, то для синкліналі "(Н)« 0, для антикліналі 
(Н)» 0. В обох випадках маємо Й(Н)- - ра «0. 
ся 


Тому ГІНУ"(«Н)» 0 для синкліналі і Г(НУ"(Н)« 0 для 
антикліналі. У обох випадках для визначення 7. і 2 
використовується процес (7). Спроби обійтися грубими 
наближеннями 2: і 2, обчисленими не за схемою (7), 


успіху не мали, хоча у ряді випадків призводили до до- 
брих результатів. 





1 Р ЗРРЕРИНИНО Ирина орат ори 


Рис. 2. Моделі контактів, описувані функціями з класу "антикліналей" (9) и "синкліналей" (10) 





2 
2 -4Н агсід- - 
Н 


Так, якщо покласти 9 -2сіпо----з 
с": 


С і 
-47 агсід-, то для синкліналі 2-Нчи, Н, 2, и»д0, 
7 


с2 а, о (знак рівності - при с-»с) і, отже, 


4Н агоід 2--А(Н чи) агсід- 7- з4На, -А(Н иа, «0. 
Н Най 


Тому 


- Т 2 мч 





сан? 9 
ІП----- 2, 
сова" 2 
9 
са; «(с НО 2, (11) 


38 
и «4(с" --Н?)е 20 - с" -Н 


Для антикліналі 2-Н-и, Ніи»д0, а,са, і, отже 


АН агсід р -АН -и) агсід -4На, - 42а, 20. Тому 
7 


етан « 9 

сова с 

При невеликих значеннях и (порівняно з Н ) можна 
вчинити інакше, а саме, взяти величини 


2Н 
піс? «(«(НаийУу1-т(с? --Н?) є -----з 
Ї ( у1-1п( ) чна 


9 
и«Н-Ч(с"--Н?)е 2 -с". (12) 


С С С СН 
Н--иу)агсій--- Н агсід -- є | агсід-- --з---- | 
ЗГУНУ Зн Зн Р 

с. 
для синкліналі. Тому 4 - -4 агеа 2О і, отже, 


г-)ау/а аг (13) 


оскільки для синкліналі )9«0. Для антикліналі, хоча 
9»0,але 2-Н-и, тому формула (13) є вірною. За- 


уважимо, що Г(Н)-9, РН) - 4 аг і тому (13) є 


першим наближенням в ітераційному процесі (7). 
Інший спосіб обчислень базується на визначенні рі- 
зних моментів кривої 2 - 2(2). У простому випадку, Згі- 


дно теореми Гауса Гасоах -2лт,де т -аномальна 
маса, зосереджена в області, обмеженій кривими 
2, - 22), 2, --Н, ас5«ф, вісь 07 спрямована вгору. 


Через це маємо 


т- (одбаєз | знає - ((аннуає о -1. 








Звідси 2. 2 -Н,2х -Н або 
-а -а 
1 о 
2 хуах-н, 
37 оцр-а) 1409 
| | 9(хуах -Н (14) 
" 2х(-а)?, 
У ЦИХ формулах всюди вважаємо, що 


о а 
Г асоах «Ге(мах. Для обчислення цього інтегралу 


-00 С 


отримано відповідний вираз у квадратурах. 
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Висновки. Нові аналітичні конструкції (6) для обчи- 
слення багатошарових контактів при їх чисельному мо- 
делюванні (7) за способом Ньютона |5| швидше збіга- 
ються порівняно із класичними ітераціями (1). Їхню стій- 
кість на даних великої розмірності доцільно перевірити 
на польових даних. 

Спроби обійтися грубими наближеннями, обчисле- 
ними не за (7), а за (11), (12), загалом, успіху не мали: 
збіжність на порядок менша, ніж у |1), математично ко- 
ректні результати мають сумнівний геологічний зміст. 
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ОМІОУуМЕЗ5 ОЕ АРРКОХІМАТІОМ САТ СЧІАТІОМ5 ЕОК МОСТИ АУЕВБЕО ОЕМ5ІТУ ІМТЕКБЕАСЕЗ 


Тпе доаїз ої пе рарег аге іо обіаїп таїпетаїіса! соп5ігисіїоп5 їог деоіодісаї! об)есіз, 5исП а5 5упсіїпе5 апа апіїсіїпе5; іо 5зиб5іапіаїе пе 
ипідиепе55 ої пе іпуег5е ргоріет у/пеп гепоуаїйпод апаїуйїса! тоавєеіз Гог пе Погігопіаїіїу Іауегеа деоіодіса! теаїа у/їП 5еуега! аеп5іу іпіеіТасе5 
іп сопіасі 5игтасез5 ргеаейпеа Бу СПогпіу; апа іо у пе іесппідие5 аеуеіореа Гог іПеїг йегайуєе саїсиіайоп. 

А сотрфіпайоп ої іпез5е іо тоаєі5 аеуеіор5 а пеум/ апа тоге ассигайе арргоаспП іо угауітеїгіс іпуег5е ргобіет5 їог іпе сопіасі іпіегТасе. ТПі5 
ресотез песе5загу іо ітргоуе 5іапаага її ргоседаийгез у/Ппеп 50Іміпд іпуег5е ргобіетз5 іп бгаміїу апа тадпеїїіс Пеіаз5. 

Тпе іпуег5е ргобіет ої їпе аеп5іїу іпіегТасе іп пе Погігопіаї!їіу Іауегесаі деоіодіса! теаїа у/іїп 5еуега! аеп5іїу іпіеіТасе5 із сопіїпеа іо пе 
50Іціоп ої пе попіїпеаг іпіедга! едиайоп паї аезсгібез пе сопіасі 5ипасе гезітістеа Бу пе дімеп сопз5іапі азутріоїез у/їїпіп пе ріапаг гедіоп. 

5, іпі5 таке5 сотриіайоп тоге сотріїсаїеа Ббесаийзе ої пе ргобіет ої едиїмаіепсу 50Ійіоп5. 

Ту/о Пе!а з5ерагайоп Шеогет5 аге ргоро5беа йог іНні5 тоає! - опе йог 5еуегаї! 1-соппесіеа уоіитез апа апоїпег опе бог пе поп-сго55еа Іауегз. 
Тпе Шеогет5 ої ипідиепез55 аге Биії оп пе іеогетз5 ої Пе/а 5ерагайоп епарбіїпд іпе 50ійіоп ої іпе іпуег5е ргобіет Бу пе зиттагу ехіета! 
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дгауіу Пе!а ої п орбіесіз (оге Боаїе5, Іауег іпіеіТасез еїс.) іпгоидП пе 50Ійіоп ої Ше іпуег5е ргоріет Гог 5ерагаїе орб)есіз - Бу Пе арргоргіаїе йїе!а 
уашез5 гот іпезбе деоіодіса! оріесіз. 

Тпе питегіса! зспетез5 їог пе аеїїпійоп ої іпе іпійа! арргохітайоп ої іпе аеп5іїу іпіеіТасе іп іпе тийіїауегеа деоіодіса! теаїа аге 5іаїеаїі. 
Тпезе аідогііпт5 Гогтаїіу соіпсіае у/йНіп пе їїг5і йегайіоп. ТПеге аге а!5о ргорозеа апаїіодісаї іесппідие5 Баз5еа ої пе СПебузПпем йегайоп 
соп5ігисіїоп ог пе йегаїіуе 5ресіїїсайоп ої пе бепаміог ої пе сопіасі а5утріоїе5. 

Тпеге у/еге тоавіеа з5упіПеїїіс іпійа! арргохітайоп5 ої 5упсіїпез5 апа апіїісіїпе5 бу Шезе адогійпт5. Ап аПйегпайуе саїсиши5 теїпоа ог ії і5 
роіпіеа оці, м/нісП із ба5еа ироп Не аеїїпійоп ої пе аїйегепі тотепіз ої іПе іпіегТасе сигуе5. Бог Пе іпіедуга! саісиіайоп іпеге і5 обіаїпесі ап 
арргоргіаїе ехрге55іоп іп пе їїпіїе диаагаіигез. 

Моавєіїпу ааїа 5Поу/ іПаї пем/у апаїуйіса! соп5ігисіоп5 Гог Ше саісиіайоп ої їпе тийауегей сопіасі іпіеіТасе5 у/іїпіп їпеіїг Мем/опіап 
питегіса! арргохітайоп сопуегдуе тоге диісКіу іп сотрагі5оп у/іїП сіа55іс іесПппідие5 їог пе сопіасі аеїїпійоп. ТПеїг іпуагіабіїйу їог Ше бід 
Чітеп5іоп Йе!а ааїа 5Поиіа ре іе5іес! оп пе геа! теазигетепі5. Мо айетріз5 іо арріу гойдп арргохітайоп5 у/еге 5иссе55їйі: сопуегдепсе у/а5 
соп5іавегарбіу Іе55 ап іп ргеуіои5 са5е5, апа, Ббезіае5, іпеге у/а5 а гаїпег атбідиоц5 деоіодісаї! таїіпіепапсе. 

Кеу у/огаз: роіепіїа! іпеогу, апаїіуйса! тоаєі, сопіасі ргобіет, сіа55е5 ої аеп5іїу іпіеіТасе5, угауйу йеій5 5ерагайїіоп, тоавєіїпд. 
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ОБ ОДНОЗНАЧНОСТИ ВЬІЧИСЛЕНИЯ ПРИБЛИЖЕНИЙ МНОГОСЛОИСТЬІХ ПЛОТНОСТНЬХ КОНТАКТОВ 


Цель статьи - получить математические конструкции для геологических обьектов типа синклиналей и антиклиналей, обосно- 
вать единственность обратной задачи восстановления аналитических моделей горизонтально-слоистой геологической средь с 
несколькими плотностньтми границами раздела для зтих конструкций в наперед определенном классе Черного контактньтх поверх- 
ностей и апробировать разработанную методику для их итерационного вьгічисления. Совокупность зтих двух моделей образует 
новую уточненную постановку обратной задачи гравиметрии для контактной поверхности. Зто необходимо для улучшения извест- 
ньїжх процедур подбора в решений обратньгх задач гравитационньтх и магнитньтх полей. 

Обратная задача определения контакта в горизонтально-слоистой среде с несколькими плотностньгми границами сведена к ре- 
шению нелинейного интегрального уравнения, которое описьваєт контакт, ограниченньтй заданньми постоянньїюми асимптотами в 
плоской области. Однако в подобной постановке практика вьгтчислений усложняется проблемой зквивалентности решений. 

Для зтой модели приведень две теоремь разделения полей - для случая нескольких односвязньтх обьемов и для случая несколь- 
ких непересекающихся слоев. Теоремьг единственности основань на теоремах разделения полей, которьгте позволяют свести реше- 
ние обратной задачи по суммарному внешнему полю п обьектов (рудньіх тел, границ раздела слоев) к решению обратной задачи для 
отдельньх обьектов - по значениям полей от зтих обьектов. 

Указаньт численньге схемьт для определения начального приближения плотностного контакта в многослоистой геологической 
среде. Зти алгоритмь! формально совпадают на первом шаге итераций. Аналогичньге схемьт на основе итерационной конструкции 
Чебьішева предложеньт и для итерационного уточнения поведения "асимптот" контакта. 

Осуществлено моделирование синтезированньтх начальньтїх приближений "антиклиналей" и "синклиналей" по зтим алгоритмам. 
Указан альтернативньтй способ вьгтчислений, которьгй базируется на определений различньгїх моментов кривой контакту. Для вьічис- 
ления интеграла получено соответствующее ввіражение в конечньтх квадратурах. 

По результатам моделирования вьгявлено, что новьгї аналитические конструкции для вьїтисления многослоистьтх контактов 
при их численном моделирований способом Ньютона скорее сходятся в сравнений с классическими методами вьїчислений контакта. 
Их устойчивость на данньїх большой размерности целесообразно проверить на полевьіх данньгх. Попьштки обойтись грубьгми при- 
ближениями успеха не имели: сходимость на порядок меньше и довольно сомнительное геологическое содержаниє. 

Ключевьг слова: теория потенциала, аналитическая модель, контактная задача, классьт: контактньгх поверхностей, разделение 
гравитационньтх полей, моделированиєе. 
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У статті з феноменологічних позицій концептуально охарактеризовано п'ять основних (або фундаментальних) прин- 
ципів, на яких базується методологія дослідження особливостей будови земної кори на підставі використання дистанцій- 
них аерокосмічних технологій. Головний (перший) принцип стосується усвідомлення процесу одночасної взаємодії (параге- 
незу) фізичних сил протилежного спрямування. Це сили притягання |і відмінні від них сили відштовхування. 

Другий принцип полягає у визначенні двох груп геологічних рушійних сил (внутрішніх і зовнішніх) у виникненні та роз- 
витку структур літосфери. Третій принцип застосування матеріалів аерокосмічних зйомок у геологічних цілях базується 
на позиціях існування в літосфері поряд з пластичними (квазіпластичними) - крихких (квазікрихких) деформацій гірських 
порід. Четвертий принцип стверджує існування активних процесів перманентно-перервного енергомасообміну в природ- 
них геосистемах планети, як проявів особливостей існування і форм руху матерії. 

П'ятий принцип констатує, що специфічні сучасні ландшафти та притаманні їм аномалії, будучи геоіндикаторами 
внутрішньої структури геологічного субстрату і пов'язаних з ним енергомасообмінних процесів, диференціюються за 
спектральними характеристиками відбитого, розсіяного, поглинутого й емітерного електромагнітного випромінювання 
в різних діапазонах хвиль. Як наслідок, вони розрізнюються за цими характеристиками, а також можуть бути відображе- 
ними і розпізнаними на даних дистанційних аерокосмічних знімань. 

Ключові слова: методологія, земна кора, геологічні рушійні сили, дистанційне зондування Землі, ландшафтні гео- 
індикатори. 


Мета та задача статті. Дистанційні знімки Землі з 
літальних космічних або повітряних апаратів є важли- 
вим джерелом даних при вивченні блокової будови її 
кристалічної основи та осадової товщі, уточненні роз- 
ташування кільцевих структур, прогнозі нафтогазо- і 
рудоперспективних об'єктів і просторового розподілу 
їхніх проявів і родовищ, дослідженні геодинаміки та 
моніторингу сучасних екзогенних процесів тощо 


1; 5; 6; 11; 12; 14 та ін... Використання даних дистан- 
ційного зондування Землі (ДЗЗ) у геологічних дослі- 
дженнях аргументується низкою їхніх позитивних якос- 
тей, важливими з них є: актуальність, оперативність 
отримання, об'єктивність у відображенні об'єктів і про- 
цесів, велика оглядовість, інтегруюча здатність (приро- 
дна генералізація), своєрідна "рентгеноскопічність", 
повторюваність, можливість отримання дистанційних 
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